home *** CD-ROM | disk | FTP | other *** search
/ Pascal Super Library / Pascal Super Library (CW International)(1997).bin / MATH / NRPAS13 / SNCNDN.DEM < prev    next >
Text File  |  1991-04-29  |  716b  |  28 lines
  1. PROGRAM d6r30(input,output,dfile);
  2. (* driver for routine SNCNDN *)
  3. VAR
  4.    cn,dn,em,emmc,sn,uu,val : real;
  5.    i,nval : integer;
  6.    txt : string[26];
  7.    dfile : text;
  8.  
  9. (*$I MODFILE.PAS *)
  10. (*$I SNCNDN.PAS *)
  11.  
  12. BEGIN
  13.    glopen(dfile,'fncval.dat');
  14.    REPEAT readln(dfile,txt) UNTIL (txt = 'Jacobian Elliptic Function');
  15.    readln(dfile,nval);
  16.    writeln(txt);
  17.    writeln('mc':4,'u':8,'actual':16,'sn':13,
  18.       'sn^2+cn^2':15,'(mc)*(sn^2)+dn^2':18);
  19.    FOR i := 1 to nval DO BEGIN
  20.       readln(dfile,em,uu,val);
  21.       emmc := 1.0-em;
  22.       sncndn(uu,emmc,sn,cn,dn);
  23.       writeln(emmc:5:2,uu:8:2,val:15:5,sn:15:5,
  24.          (sn*sn+cn*cn):12:5,(em*sn*sn+dn*dn):14:5)
  25.    END;
  26.    close(dfile)
  27. END.
  28.